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ABSTRACT 

The most luminous GRBs can be explained in terms of models involving 
stellar mass central engines only if the ejecta are beamed. As was pointed 
out by Rhoads (1997), the dynamics of the blast wave, formed by the beamed 
ejecta sweeping the external gas, can be significantly modified by the sideways 
expansion. This is because in this case the surface of the blast wave increases 
faster than just due to the radial divergence and so the blast wave deceleration 
rate increases faster. According to analytical estimates, the effect becomes 
important shortly after the bulk Lorentz factor of the blast wave drops below 
the inverse of the initial opening angle of the beamed ejecta and is accompanied 
by a sharp break in the afterglow light curve. 

However, our numerical studies, which follow the dynamical evolution of 
the blast wave, the evolution of the electron energy distribution, and take into 
account the light travel effects related to the lateral size of the source, show 
that the break of the light curve is weaker and much smoother than the one 
analytically predicted. A prominent break emerges only for a model without 
sideways expansion. 

Subject headings: gamma rays: bursts 
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1. INTRODUCTION 

Beaming of relativistic ejecta in GRBs has been postulated by many authors in order 
to ease the GRB energy budget (see, e.g., Meszaros, Rees, & Wijers 1998 and refs. therein). 
There are basically two ways to verify the beaming observationally: one is statistical 
and is based on counting the afterglow like transient sources and comparing their rate 
with the GRB rate, and the second one is related to the beaming effects predicted to be 
imprinted in the afterglow light curves of individual objects (Rhoads 1997). Applying 
the first method to the X-ray transient sources, Grindlay (1999) found that results are 
consistent with no beaming differentiating the GRB and X-transient rates. However, as 
was pointed out by Woods & Loeb (1999), the conclusive results about excess (or its lack) 
of X-ray transients over GRBs must wait for much more sensitive future instruments. This 
is because statistically significant contribution to the excess of X-ray transients over GRBs 
is expected to be provided only by weak X-ray transients, those representing afterglows 
phases when the bulk Lorentz factor of the radiating shell drops below the inverse of its 
angular size. Similar studies can be performed also in optical and radio band (Rhoads 1997; 
Woods & Loeb 1999). 

In individual objects, the beaming related effects are expected to be imprinted in the 
optical and X-ray afterglow light-curves. The lateral expansion of the shocked, relativistic 
plasma causes that at some moment the front of the blast wave starts to increase faster 
than due to the cone-outflow (Rhoads 1997). Due to this the blast wave begins to decelerate 
faster than without the sideways outflows and this produces a break in the light curve, the 
sooner the larger the beaming factor is. Such a break is claimed to be present in the light 
curve of GRB 990123, the most energetic GRB up to date (Kulkarni et al. 1999). Sari, 
Piran, & Halpern (1999) speculate that afterglows with very steep light curves are highly 
beamed. Possibly the break in such objects is not recorded because it took place before the 
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optical follow-ups. 

As now, all theoretical studies of the light-curve breaks are analytical and are based 
on: power-law approximation of the blast wave dynamics, broken power-law approximation 
of radiation spectra and on "on-axis" relation between the observed flux and the emitted 
flux (Rhoads 1977, 1999; Kulkarni et al. 1999; Sari, Piran, & Halpern 1999). In this paper, 
we treat the dynamics using the prescription given by Blandford and McKee (1976). The 
evolution of radiation spectrum is calculated exactly, by computing the time evolution 
of electrons from continuity equation and by computing the observed luminosity through 
integrating the emitted radiation over the "t = const" surfaces. Our results show that the 
change of the light curve slope is significant, but smaller than predicted analytically. And, 
what is more important, the light curves steepen very slowly, so that it is very difficult 
to talk about the specific time location of the break. In order to demonstrate better the 
beaming effect, we compare our results with the spherical case. We also show how the light 
curve should look like if there is no lateral expansion. 

In §2 we collect equations, which are used to compute the blast wave speed, evolution 
of electrons, and afterglow light-curves. In §3 we present results of our numerical studies of 
afterglows produced by beamed ejecta, and, in §4 we compare them with simple analytical 
estimations. 



2. BASIC EQUATIONS 
2.1. Dynamics 

The deceleration of a blast wave is described by the following equations (Blandford & 
McKee 1976; Chiang & Dermer 1998): 
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dM dm r _ 

■^ = ^[r- erajee( r-i)], (2) 

and 

dm/dr = Qjr 2 p = 2nr 2 (l — cos9j)p, (3) 

where T is the bulk Lorentz factor of the blast wave, M is the total mass including internal 
energy, r is the distance from the central engine to the blast wave, dm is the rest mass 
swept up in the distance dr, p is the mass density of the external medium, e e is the fraction 
of dissipated energy converted to relativistic electrons, e ra d is the fraction of electron energy 
which is radiated, and 9j is the angular size of the blast wave. This angular size is not 
constant but increases due to thermal expansion (Rhoads 1997), 

e j = ± = o j0 + ± (4) 

where the speed of the lateral expansion, v[, is assumed by Rhoads (1999) to be equal to the 
sound speed in the relativistic plasma, c s = c/v^3, but considered by Sari et al. (1999) to be 
relativistic. Noting, that the plasma in the blast wave is continuously loaded by the fresh 
gas which initially doesn't have any lateral bulk speed, one can expect that in reality v[ 
does not reach relativistic value and sets up somewhere between c s and flrc, and in general 
depends on r, and on 9j. 

2.2. Electron energy distribution 

We assume that the electrons are injected with the power low energy distribution 

Q = Kr v , (5) 
with the minimum energy of injected electrons 

= e . ( r - ik v_^2 (6) 

m e p-1 
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The maximum energy of injected electrons for a given magnetic field, B', is assumed to be 
given by (de Jager et al. 1996): 

f B ,\ -1/2 

7.^4X10'- (7) 
Normalization of the injection function, K, is provided by 



Le,inj = / Q~fm e c 2 d-i = e e —f^ , 

J"fm 01 



where 



dE' acc dr dE' acc dr dm 2 2 3 

-IT = It— = dv^ {T - l) = n ^"^ T - 1 ) c • < 9 > 



is the rate of accreted kinetic energy, dr = cfirFdt', fir = Vr 2 — 1/r, and t' is the time 
measured in the blast wave comoving frame. 



The evolution of the electron energy distribution is given by the continuity equation 

-H-(n,P)+q, (10) 



ON, _ d AT d 1 
dr d'j \ 7 dr 



where 



^ = -/(r)7 2 -^, (11) 
dr r 



are the electron energy losses. In both equations above, the derivatives over comoving 
time t' have been replaced by the derivatives over the distance r, according to the relation 
d/dt' = c/3rTd/dr . The first term on the rhs of Eq. (|ll|) represents synchrotron plus 
Compton energy losses, i.e., 

^=£^S (i+<M ' (12) 

where u' B = B' 2 /8 it is the magnetic energy density, and u' s is the energy density of the 
synchrotron radiation, both as measured in the blast wave frame. The second term on the 



rhs of Eq. (11) represents the adiabatic losses. The parameter g depends on the geometry 
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of the expansion; for 2-dimensional (lateral) expansion g = 2/3, and for 3-dimensional 
expansion g — 1. 

We calculate the magnetic field following Chiang & Dermer (1999) 

B' 2 

u' B = — = e B K,pc 2 T 2 , (13) 

where k is the compression ratio and es parameterizes the departure of the magnetic field 
intensity from its equipartition value. 



2.3. Synchrotron spectrum 

The evolution of the synchrotron spectrum in the blast wave frame is given by 

L' syn y(r) = /iV 7 (r)P(,' )7 )d 7) (14) 

where P(v', 7) is the power spectrum of synchrotron radiation of a single electron in 
isotropic magnetic field (see, e.g., Chiaberge & Ghisellini 1999). 

The apparent monochromatic synchrotron luminosity as a function of time (a light 
curve) is calculated from 

LsynAt^obs) = JJ^ ~ — dcos6d(j), (15) 

where T> = 1/T(1 — /^pcos^) is the Doppler factor of the blast wave at the angle 0. The 
coordinates (9, <p) are chosen so that the observer is located at 9 = {9 = 9 Q b s ) and the jet 
axis is at 9 = 9 obs . The integral is taken over the surfaces 

f (1 - B T cos 9) , 

t = / ^ '- dr = const, (16) 

J cpr 

enclosed within the blast wave boundaries, Qj. 
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2.4. Inverse-Compton radiation 

We assume hereafter that cooling of relativistic electrons is dominated by synchrotron 
radiation, i.e. that u' s <C u' B . This condition will be verified and discussed in Appendix A. 

3. RESULTS 

We have used the following model parameters of the afterglow model in our calculations: 
initial energy per solid angle, E /£lj = 10 54 ergss _1 /47r; T = 300; 9j = 0.2; k = 4; 
p = m p /lcm 3 ; e e = 0.1 (the quasi-adiabatic case); e B = 0.03; p = 2.4. The parameters 
were not chosen to fit any specific observations, but rather to demonstrate the difference 
between simple analytical predictions and self-consistent numerical calculations regarding 
the beaming effects in a light-curves. 

In Fig. 1 we present the dependence of a bulk Lorentz factor of the blast wave on its 
distance from the central engine. We show three solutions for three different values of v[: 
v[ = (thin line); v[ = c/\/3 (solid line); and v[ = c (dotted line). For v[ = (— > 9j = const) 
and r <r< r nr , the bulk Lorentz factor is well approximated by 

r~r (r /r) 3/2 , (17) 

where 

is the radius where deceleration of the GRB ejecta by sweeping of interstellar gas starts to 
be efficient, and 

~3.4xl0 18 cm, (19) 

is the radius above which the blast wave becomes nonrelativistic. Fig. 1 demonstrates that 
steepening of the T(r) curves due to lateral outflow is very smooth, without any sharp 
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break like the one predicted analytically to take place at a distance, at which T drops below 
l/Ojo, i.e. at 

r D ~ r (r % ) 2/3 - 1.8 x 10 18 cm . (20) 

In Fig. 2 we show the radial dependence of the rate of kinetic energy accreted by the 
blast wave, dE' acc /df (see Eq. ||). As one could expect, the larger the lateral outflow speed, 
the larger the accretion rate is. The steepening of curves at large r is due to transition from 
the relativistic regime (r > 2) to the nonrelativistic regime, where dE' acc /dt' is significantly 
reduced, and becomes oc (r — 1) (see Eq. |9|). When divided by e e , curves in Fig. 2 illustrate 
also the r dependence of injection luminosity of relativistic electrons (see Eq. ||). 

In Fig. 3 we show time evolution of the electron energy distribution, N 7 , multiplied by 
7 2 . The curves are calculated at such values of the radius r, from which the signal produced 
on the axis 6 = is reached by the observer t = 1,10, 10 2 , 10 7 seconds after "the signal" 
from r = 0. The relation between r and t is 

t-rZzM^fi*. (21) 

Jo cpr Jo 2cT z 

The peak positions of the iV 7 7 2 curves mark Lorentz factor of those electrons which carry 
most of leptonic energy at a given distance. For injection spectral index 2 < p < 3, the peak 
is located at 7 m given by Eq. (|5[), and this is the case in our model. Another characteristic 
energy is 

7c = 6.1xl0 20 ^4, (22) 
EgKp rl 

below which the time scale of electron energy losses due to synchrotron radiation is longer 
than the dynamical time scale. We present the dependence of 7 m and 7c on the radius r in 
Fig. 4. 

We can see from Fig. 4, that for the first 5 curves presented in the Fig. 3, 7m > 7c . 
In this case, in accordance with the analytical predictions, the electron spectra at 7 > 7 m 
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are well described by the power-law function, iV 7 oc 7 _s , with the index s = p + 1. For 
7 C < 7 < 7 m , analytical crude estimations predict s — 2, which in our plot should be 
represented by horizontal lines. This, however, is expected to be true only for 7 C <C 7 < j m . 
In our model the ratio jm/lc is not large enough to provide space for s = 2 and there is a 
smooth transition to very hard low energy tail reached by electrons due to adiabatic losses. 
For 7 C > 7 m , which is the case for the top 4 curves in the Fig. 3, the predicted electron 
spectra should have a slope s = p + 1 for 7 > j c , and s = p for 7 m < 7 < j c . The former is 
seen, but the latter, again, due to narrow range between 7 C and 7 m doesn't apply. Instead, 
there the log-energy distribution is curved, smoothly joining the high energy portion of 
the electron spectrum with its low energy adiabatic part. Let us note, that very steep low 
energy tails of the curves on top of the plot result from the fact that there is not enough 
time for electrons to drift adiabatically to lower energies. Note also, that the details of the 
low energy parts of the electron energy distribution are not important, because contribution 
of electrons from these parts to the observed radiation is negligible. 

In Fig. 5 we present the observed radiation spectra computed for the same sequence 
of t, as electron energy distributions shown in Fig. 3. It should be pointed out, however, 
that unlike in simple analytical calculations, they are computed by integration of electron 
radiation from t = const surfaces (see Eq. [IB]), i.e. taking into account light travel differences 
between photons emitted at different #'s. The observed radiation spectra are peaked around 
hv ~ T^B / B cr m e c 2 , where rhs quantities are calculated for r given by Eq. ([16]) and 
B cr = 2ixm 2 e c' 1 'he ~ 4.4 x 10 13 Gauss. As we can see from Fig. 5, the high energy and the 
low energy parts of the observed radiation spectra are well described by power-law functions 
L v oc v~ a , with a = p/2 = 1.1 and a = —1/3, respectively. The former is produced, as 
predicted analytically, by electrons with 7 > Max[7 m ; 7J, the latter represents the low 
energy synchrotron radiation of electrons with energies 7 < Min[7 c ; 7 m ]. These high and 
low energy spectrum portions are joined very smoothly without showing any intermediate 
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piece of the power-law spectrum. This smoothing of the observed radiation spectra results 
mainly from the fact that the observed radiation at any given moment is contributed by 
radiation from t = const surfaces, i.e. from different radii. 

The light curves, computed for v = 4.2 x 10 14 Hz and v = 2.5 x 10 17 Hz, are shown in 
Fig. 6. In this calculation we used v[ = cj and considered two different locations of the 
observer, 9 = and 9 = 0.28. The latter case is for the observer located outside the initial 
ejecta cone. 

In order to demonstrate better the beaming effect and its dependence on v' h we plot 
in Fig. 7 four optical light curves; three for different values of v[: 0, c/VS, and c, and the 
fourth one for the spherical outburst. We can see, that for models with lateral expansion 
the steepening of the light curves is extended over more than two time decades and, down to 
nonrelativistic regime, doesn't reach analytically predicted slope (5 = p [L v oc t^ 13 ). Sharp 
break is found only for v[ = model, and, it emerges shortly after t(r D ), as theoretically 
predicted. 

We should note here, that our calculations of the model with the outflows expansion 
are not fully consistent, because the Doppler factor includes only the radial component of 
the bulk motion. This, however, is expected to affect only the results of the vi = c model, 
where we overestimate the radiation contributed from the blast wave edge. 

4. DISCUSSION AND CONCLUSIONS 

The afterglows provide exceptional opportunity to study whether and how much the 
GRB ejecta are beamed. As predicted by Rhoads (1997), the beamed outflows should 
diverge from the cone geometry while decelerated by sweeping up the external gas. The 
sideways outflow of the shocked relativistic plasma increases the front of the blast wave 
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leading to a faster deceleration. Rhoads (1997) showed, using simple analytical analyses, 
that this should be imprinted in the light curve as a break around t(r D ), i.e. when T 
drops below (q/c)/6>. There the light curve should steepen, changing the slope from 
(3 = (3p - 2)/4 (for 7 > Max[ 7c ; 7m ]) or from (3 = 3(p - l)/4 (for 7m < 7 < 7c ) (Sari, Piran, 
& Narayan 1998) to /3 = p (Rhoads 1999). Our numerical results agree only qualitatively 
with these predictions. The steepening does occur, however, the slope change is smaller (to 
about 2.0 instead of p — 2.4) and is extended over more than two decades of the observed 
time. 

There are two reasons why, contrary to simple analytical estimations, the distinct break 
does not emerge in our calculations. First, as is shown in Fig. 1, the dynamics of the blast 
wave is affected by the lateral outflow very smoothly over the whole deceleration phase, and 
not just around r D (note, that at r D the blast wave area is already almost 4 times larger 
than it would be without the sideways expansion). Second, the observed radiation at any 
given moment t is contributed by the plasma which at larger 9 emits radiation from smaller 
r, and noting, that at smaller r plasma is moving faster and radiating stronger than at 
larger r, the contribution of the off- axis plasma to the observed radiation is larger than in 
the case of radiation contribution taken from r = const surfaces as calculated analytically. 

The steepening of the light curve is predicted also for the beamed ejecta without the 
lateral outflows. In this case, because no change of dynamics and small light travel effects 
at r > T} } (note that there the Doppler cone becomes narrower than the ejecta cone), the 
break is very well located, just around the time the Y drops below 1/9 j, and the light curve 
steepens by A(3 = 3/4, in accordance with analytical predictions (see, e.g., Meszaros & 
Rees 1999). 

It should be emphasized, however, that our treatment of the dynamics with the 
sideways expansion is based on the approximation, that at any r the material is uniformly 
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distributed across the blast wave. In reality, the lateral outflow can create 9 dependent 
structure, with the density of the swept material and the radial bulk Lorentz factor 
decreasing sideways, and in this case the break in the light curve may become more 
prominent. 2D- hydro relativistic simulations are required to verify this. 
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A. INVERSE COMPTON COOLING 



The ratio of inverse Compton luminosity to synchrotron luminosity is given by 

if = t> (A1) 



where 



This manuscript was prepared with the AAS L^TjrjX macros v4.0. 
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L' s = (l- Vc )e rad e e -^, (A3) 
dE' acc /dt' and u' B are given by Eq. (g) and (0), respectively, and r\ c = L' C /(L' C + L' s ). Using 



all these relations in Eq. ( |A1| ), we find that for r 3> 1 

L' c (1 - r] C )e rad e £ 



L' s 2ne B 
and noting that L' c /L' s = r] C /(l — r] C ), we obtain 



(A4) 



where 



and for x C 1, ?/c - X' 



(1 + 2 X ) - 

= ^ , (A5) 



X = e ~^- (A6) 



For 7 m > 7 C practically all energy converted to electrons is radiated (e ra d — 1) and 
the inverse Compton is energetically negligible (jj C -C 1/2) if e B > 10~ 2 . For 7 m < j c , the 
luminosity peaks at v c and 

e ra ^(^) P2 . (A7) 
In the latter case, the inverse Compton is energetically not important, if 

^{^-f. (A8) 

7m \ e B J 

One can easily check, using the above criteria and Fig. 4, that for our specific model the 
inverse Compton process does not dominate electron cooling at any moment. 

It should be noted, however, that the inverse Compton process can be imprinted in the 
afterglow light curves, even if the Compton cooling is less efficient than the synchrotron 
one. At the moment when the Compton component drifts down to the observed band, the 
light curve is expected to flatten. Chiang & Dermer (1999) demonstrated that this effect 
can be visible in the X-ray light curves; in the optics it appears very late and is already too 
weak to be observed, especially if overshined by the host galaxy. 
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FIGURE CAPTIONS 

Fig. 1. — The bulk Lorentz factor of the blast wave as a function of a distance from the 
central engine. Thin line: for 9j = const = 9jq; solid line: for 9j = 9jo + l/y/3T; dotted line: 
for 9j = 9jo + 1/r. The numbers along the curves show the observed time as measured in 
seconds and given by Eq. (^Tj). 

Fig. 2. — The accretion rate as a function of a distance from the center, calculated for the 
same models as Fig. 1. 

Fig. 3. — Evolution of the energy distribution of relativistic electrons. From the bottom to 
the top, the curves are for t: 1, 10, 10 7 seconds. 

Fig. 4. — The minimum injection electron energy, j m , and the "cooling" electron energy, j c , 
as a function of time. 

Fig. 5. — The evolution of the apparent radiation spectra. The numbers mark the observed 
times. 

Fig. 6. — The afterglow light curves: the dotted lines are the X-ray light curves and the solid 
lines are the optical light curves. The thick lines are for 9 Q b s = while the thin lines are for 
Gobs = 0.28. 

Fig. 7. — The light curves in the time range, where the beaming effects are strongest. The 
thin solid line is for 9j = 9jo; the thick solid line is for 9j = 9jo + l/^/3T; the dotted line is 
for 9j = 9jo + 1/r; and the dashed line is for 9j = n (spherical case). 
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